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Abstract 

We analyse the renormalizability of the sine-Gordon model by the example of 
the two-point causal Green function up to second order in a r (M 2 ), the dimensional 
coupling constant defined at the normalization scale M, and to all orders in f3 2 , the 
dimensionless coupling constant. We show that all divergences can be removed by 
the renormalization of the dimensional coupling constant using the renormalization 
constant Z\, calculated in (J. Phys. A 36, 7839 (2003)) within the path-integral 
approach. We show that after renormalization of the two-point Green function to 
first order in a r (M 2 ) and to all orders in j3 2 all higher order corrections in a r {M 2 ) 
and arbitrary orders in /3 2 can be expressed in terms of Qf p h, the physical dimensional 
coupling constant independent on the normalization scale M . We calculate the Gell- 
Mann-Low function and analyse the dependence of the two-point Green function 
on a p h and the running coupling constant within the Callan-Symanzik equation. 
We analyse the renormalizability of Gaussian fluctuations around a soliton solution. 
We show that Gaussian fluctuations around a soliton solution are renormalized like 
quantum fluctuations around the trivial vacuum to first orders in a r (M 2 ) and (3 2 
and do not introduce any singularity to the sine-Gordon model at /3 2 = 8ir. The 
finite correction to the soliton mass, coinciding with that calculated by Dashen et al. 
(Phys. Rev. D 10, 4130 (1974)), appears in our approach to second order in a p h and 
to first order in f3 2 . This is a perturbative correction, which provides no singularity 
for the sine-Gordon model at j3 2 = 8ir. We calculate the correction to the soliton 
mass, caused by Gaussian fluctuations around a soliton, within the discretization 
procedure for various boundary conditions and find complete agreement with our 
result, obtained in continuous space-time. 
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1 Introduction 



C(x) = - d^{x)d^{x) + -^-^ (cos/?tf(x) - 1), (l.r 



C{x) = -d^{x)d^{x) + rK 2 (cos/3#{x) - 1) + (Zi - 1) ^ ' (cos/M(x) - 1) = 
~ d^{x)d^{x) + Z x (cos/3tf(x) - 1), (1.2) 



The sine-Gordon model we describe by the Lagrangian PJ E] 

±d^(x)d^(x) + ^f" 

where the field $(x) and the coupling constant (3 are unrenormalizable, ao(A 2 ) is a di- 
mensional bare coupling constant and A is an ultra-violet cut-off. As has been shown in 
j2] the coupling constant ao(A 2 ) is multiplicatively renormalizable and the renormalized 
Lagrangian reads 

~d^(x)d»V(x) + ^^(cos/M(x) - 1) + - l)^f ' 

1 a (M 2 ^ 

- d^{x)d^{x) + Z 1 ^- 2 

where Z\ = Zi(a r (M 2 ), (3 2 , M 2 ; A 2 ) is the renormalization constant [2] _ |S] depending 
on the normalization scale M. The renormalization constant relates the renormalized 
coupling constant a r (M 2 ), depending on the normalization scale M, to the bare coupling 
constant ct (A 2 ) 0: 

a r (M 2 ) = Z^(a r (M 2 ),(3 2 } M 2 ; A 2 )a (A 2 ). (1.3) 

As has been found in [2] the renormalization constant Zi(a r (M 2 ), (3 2 , M 2 ; A) is equal to 

/ A 2 \ P 2 /^ 

Z l (a r (M 2 ),(3 2 ,M 2 ;A 2 )= (— ) . (1.4) 

This result is valid to all orders of perturbation theory developed relative to the coupling 
constant [3 2 and a (A 2 ) [2J. Since the normalization constant does not depend on a r (M 2 ), 
we write below Z x = Z^/3 2 , M 2 ; A 2 ). 

For the analysis of the renormalizability of the sine-Gordon model with respect to 
quantum fluctuations around the trivial vacuum we expand the Lagrangian (jl.2|) in powers 
of i9(x). This gives 

= ^ld^(x)d^(x)-a r (M 2 )$ 2 (x)]+£ int (x), (1.5) 
where C- m t(x) describes the self-interactions of the sine-Gordon field 



C int (x) = a r {M 2 )Y,\^P 2{n - l) $ 2n (x) 



71=2 



00 (—l) n 

[Zi - 1) a r (M 2 ) J2 V-^ 20 ^ 1 ^)- 

n=l [ - AH) - 



(1.6) 



It is seen that the coupling constant a r (M 2 ) has the meaning of a squared mass of free 
quanta of the sine-Gordon field $(x). The causal two-point Green function of free sine- 
Gordon quanta with mass a r (M 2 ) is defined by 

-.A F (x;a r (M 2 )) = (0|T(^(0))|0)= / -i* * - - . (1.7) 



(2it) 2 i a r (M 2 )-k 2 -iO' 
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At x = the Green function — i A^(0; a r (M 2 )) is equal to [2] 

-iA F (0;a r (M 2 )) = — in 

where A is a cut-off in Euclidean 2-dimensional momentum space 

The paper is organized as follows. In Section 2 we analyse the renormalizability of 
the sine-Gordon model by means of power counting. In Sections 3 we investigate the 
renormalizability of the two-point Green function of the sine-Gordon field. Quantum 
fluctuations are calculated relative to the trivial vacuum up to second order in a r (M 2 ) and 
to all orders in (3 2 . We show that after renormalization of the two-point Green function 
to first order in a r (M 2 ) and to all orders in (3 2 all higher order corrections in a r (M 2 ) 
and all orders in /3 2 can be expressed in terms of a p h, the physical dimensional coupling 
constant independent on the normalization scale M. We derive the effective Lagrangian 
of the sine-Gordon model, taking into account quantum fluctuations to second order in 
a r (M 2 ) and to all orders in (3 2 . We show that the correction of second order in a p ^ 
and first order in /3 2 , i.e. 0(a 2 h j3 2 ), reproduces the finite contribution — y/aph/n to the 
soliton mass coinciding with that calculated by Dashen et al. [HI Ej. In our approach 
this is a perturbative correction, which does not lead to a singularity of the sine-Gordon 
model at j3 2 = 8tt. This confirms the absence of a singularity of the renormalized sine- 
Gordon model at f3 2 = 8tt conjectured by Zamolodchikov and Zamolodchikov jS] and 
proved in j2]. In Section 4 we analyse the renormalizability of the sine-Gordon model 
within the Renormalization Group approach. We use the Callan-Symanzik equation 
for the derivation of the total two-point Green function of the sine-Gordon field in the 
momentum representation. We show that the two-point Green function depends on the 
running coupling constant a r (p 2 ) = aph(p 2 / a ph)^ ^ where (3 2 = f3 2 /{l + (3 2 /8n) < 1 
for all (3 2 . In Section 5 we investigate the renormalizability of the sine-Gordon model 
with respect to Gaussian fluctuations around a soliton solution. We show that Gaussian 
fluctuations around a soliton solution lead to the same renormalized Lagrangian of the 
sine-Gordon model as quantum fluctuations around the trivial vacuum taken into account 
to first order in a r (M 2 ) and f3 2 . In Section 6 we discuss the correction to the soliton mass 
induced by quantum fluctuations. We show that Gaussian fluctuations around a soliton 
solution reproduce the same correction as the quantum fluctuations around the trivial 
vacuum, calculated to first orders in a r (M 2 ) and (3 2 . This correction does not contain 
the non-perturbative finite quantum correction obtained by Dashen et al. [Sid- I n our 
analysis of the sine-Gordon model such a finite correction appears only to second order 
in a ph and to first order in (3 2 (see Section 3). In Section 7 we discuss the calculation 
of the correction to the soliton mass AM S , induced by Gaussian fluctuations, within a 
discretization procedure for various boundary conditions. We show that the result of the 
calculation of AM S does not depend on the boundary conditions and agrees fully with 
that obtained in continuous space-time. In the Conclusion we summarize the obtained 
results and discuss them. In the Appendix we adduce the solutions of the differential 
equation related to the calculation of Gaussian fluctuations around a soliton. 



A 
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2 Power counting and renormalization of the sine- 
Gordon model 



As usual the general analysis of renormalizability of a quantum field theory is carried 
out in the form of power counting, the concept of the superficial degree of divergence of 
momentum integrals based on dimensional considerations [B] [E]. 

The analysis of the convergence of a given Feynman diagram G within power counting 
is done by scaling all internal momenta with a common factor A, ki — > \ke, and looking 
at the behaviour Iq ~ A w( - G ) at A — ► oo. Since we deal with a quantum field theory of a 
(pseudo) scalar field, the propagator of such a field behaves as A -2 at A — > oo. 

Let a given Feynman diagram G contain L independent loops, / internal boson lines 
and Vi n vertices with 2n lines. Since we have no vertices with derivatives of the sine- 
Gordon field, the superficial degree of divergence ui{G) of a diagram G is 

oo(G)=2L-2I. (2.1) 

The number of independent loops L is defined by [3 -[5j 

L = I + 1 - V ^ ( 2 - 2 ) 

{n} 

where the sum extends over all vertices defining the Feynman diagram G. Substituting 
(|2.2|) into (j2.1|) we can express the superficial degree of divergence w(G) in terms of vertices 
only 

cj(G) = 2-2^V 2 „. (2.3) 

{«} 

This testifies the complete renormalizability of the sine-Gordon model. 

Indeed, the vacuum energy density is quadratically divergent, since it corresponds to 
the "Feynman diagram without vertices" . Such a quadratic ultra-violet divergence of the 
vacuum energy density of the sine-Gordon model has been recently shown in [2|. Such a 
quadratic divergence of the vacuum energy density can be removed by normal-ordering 
the operator of the Hamilton density of the sine-Gordon model j2] . 

All Feynman diagrams with one vertex diverge logarithmically, and any Feynman 
diagram with more than one vertex converges. As has been shown in 2] by using the 
path-integral approach, all logarithmic divergences can be removed by the renormalization 
of the dimensional coupling constant «o(A 2 ). 



3 Renormalization of causal two— point Green func- 
tion 

The causal two-point Green function of the sine-Gordon field is defined by 

-<A(x) = i * i*Z[J] J=0 , (3.1) 
i dJ(x) % oJ(0) 
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where Z[J] is a generating functional of Green functions 
Z[J] = fm exp { i I d 2 y [C(y) + %) J(y)} 



V'd exp | i J d 2 y 



{o lll )(!))O" l ){ !l )-^ y (M 2 )0 2 { ! ,)) 4 £ ill( (//) + ^(//).7(//) j>. (3.2) 



normalized by Z[0] = 1, is the external source of the sine-Gordon field $(x). 
Substituting f)3.2j) into (j3.1j) we get 

-iA(x) = J 2W(x)tf(0) exp {i J d 2 y £ int (y)} 

x exp[ % -j d 2 y[d^{y)d^{y)-a r {M 2 )d 2 {y)]), (3.3) 

where £i n t(y) is given by (jl.6j) . The r.h.s. of ()3.3|) can be rewritten in the form of a 
vacuum expectation value of a time-ordered product 



-zA(x) = (0\T(^(x)m exp{i J d 2 y :£ int (</):})|0) ( 



(3.4) 



where the index c means the connected part, : . . . : denotes normal ordering, -$(x) is 
the free sine-Gordon field operator with mass a r (M 2 ) and the causal two-point Green 
function A F (x,a r (M 2 )) defined by (|T7j) . 

In the momentum representation the two-point Green function (J3.4j) reads 

-Uto = -i f 

d 2 Xe +'P- x {0\T(t)(x)t>(0) expji: / d 2 y :£,„,(»):)) |0> c . (3.5) 



For the analysis of the renormalizability of the sine-Gordon model we propose to calculate 
the corrections to the two-point Green function ()3.4|) (or to ()3.5|) ). induced by quantum 
fluctuations around the trivial vacuum. Expanding the r.h.s. of Eq. (j3.4J) in powers of 
a r (M 2 ) and ft 2 we determine 



-iA(x) = J^H) A (m) (x, « r (M 2 )), (3.6) 

m=0 

where (— i) A ( - m \x,a r (M 2 )) is defined by 

'TTL {* 

-iAH(x, ar (M 2 )) = — / TT Afc(0|T(^(x)^(0) : C int (y k ) :)|0) c . (3.7) 

ml J - LJ - 

J k=i 

The Green function (— i) A^(x, a r (M 2 )) coincides with the Green function (jl.7j) of the 
free sine-Gordon field. 

In the momentum representation the correction to the two-point Green function 
(— i) A^ m \x, a r (M 2 )) can be written as 



-tA^ m \p,a r {M 2 )) = J d 2 xe +t P- x {-i) A (m) (x, a r (M 2 )) = 

= % —\ d 2 xe +t P- x / l[d 2 y k (0\T($(x)i}(0) : C int (y k ) :)|0) c . (3.1 



fc=i 
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The momentum representation is more convenient for the perturbative analysis of the 
renormalization of the two-point Green function of the sine-Gordon field. 

3.1 Two— point Green function to first order in a r (M 2 ) and to all 
orders in f3 2 

The correction to the two-point Green function to first order in a r (M 2 ) and to all orders 
in (3 2 is defined by 



-zA«(x,a r (M 2 )) = i J <P yi (0\T(#(x)#(0) : CUvx) :)|0> c = 
= iar (M 2 )^ { j^f3 2 ^ J d 2 y <0|T(tf(x)tf(0) : 2n (y) :)|0> t 



Making all contractions we arrive at the expression 



aO0(O):0*(y):)|O> c . (3.9) 



00 (— 1 l n 

-i A (1) (x, a r {M 2 )) =ia r {M 2 ) V i-^ 2 ^) [2ti (2n - l)!!][-iA F (0, a r (M 2 ))] n_1 

ri ( - //)! 

x J d 2 y {-iA F (x - y,a r (M 2 ))]{-iA F (-y,a r (M 2 ))} 

+ia r (M 2 ) {Z 1 - l)J2 [ -^P 2{n ~ 1] [2n(2n - !)!!][— iA F (0, a r (M 2 ))] n_1 



x J d 2 y [—iAp(x - y, a r (M 2 ))] [—iA F (-y, a r (M 2 ))} (3.10) 
The sums over n are equal to 

£ l^Y [2n {2n ~ m [ ^ H) M °' a ^ M ^ 1 = 1 " ex P { ^ 2i M°» «,(m 2 ))}, 

n=2 ^ 

£ MT^ (n_1) [2n (2n " HM °' " r (M ' ))r 1 = " exp { \ /^ A Ho,MM 2 ))}. 

n=l \ '' ~" 

(3.11) 

Substituting (J3TT]t into (jHlDD we get 

-«A (1) (i,« r (M 2 )) = ia r (M 2 ) 1 -exp|^/3 2 zA F (0,a r (M 2 ))} 

x y d 2 yHA F (a;-y,a r (M 2 ))]HA F (-y,a r (M 2 ))] 
- ia r {M 2 ) (Z x - 1) exp{i/3 2 iA F (0,a r (M 2 )) 
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x / [-iA F (i - !/, a r (M 2 ))][-iA F (-|/, a r (M 2 ))]. (3.12) 



The r.h.s. of (|3.12|) can be transcribed into the form 

-!A (1) (i,« r (M 2 )) = ia r {M 2 ) l-Z x exp{^/3 2 zA F (0,a r (M 2 )) 



x J d 2 y [—iA F (x - y, a r {M 2 ))}[-iA F {-y, a r (M 2 ))}. (3.13) 

Using the normalization constant Z\, given by ()1.4|) . and the definition ()1.8|) of the two- 
point Green function we remove the cut-off A 



jA (1) (i,« r (M 2 )) = ia r {M 2 



a r (M 2 )\P 2 l^ 
M 2 ) 

x J d 2 y [-iA F (x - y, a r (M 2 ))][-tA F (-y, a r (M 2 ))]. (3.14) 

Thus, the renormalized causal two-point Green function of the sine-Gordon field, defined 
to first order in a r (M 2 ) and to all orders in /3 2 is given by 



- iA(x) 



iA F (x, a r (M 2 )) + i a r (M 2 ) 



a r (M 2 )\/3 2 /s- 
M 2 



x / d 2 y[-iA F {x-y,a r {M 2 ))][-iA F {-y,a r {M 2 ))}. (3.15) 



In the momentum representation the two-point Green function (|3.15|) reads 

-a r (M 2 )\/3 2 /8-i (_j) 



i A(p) 



-i 



+ 



a r (M 2 )-p 2 a r (M 2 )-p 



ia r (M 2 



M 2 



J a r {M 2 )-p 2 ' 

(3.16) 



The second term defines the correction to the mass of the sine-Gordon field 

'a r (M 2 )\^/8Ti 



5a r (M 2 ) = -a r (M 2 ) 



1 



M 2 



(3.17) 



Thus, the two-point Green function, calculated to first order in a r (M 2 ) and to all orders 
in /3 2 is equal to 



i A(p) 



a r (M 2 ) + 5a r (M 2 ) - p 2 a ph - p 2 



where a p h is determined by 



a ph = a r (M 2 ) + 6a r (M 2 ) = a r (M 2 ) l a r{M 2 )\P/** 



This gives also a r (M 2 ) in term of M and a 



P h- 



a r (M 2 ) = aJ— ) , /5 2 

\a ph / 



M 2 



/? 2 



2 ' 



1 + 



L 

8tt 



(3.18) 



(3.19) 



(3.20) 



The Green function (j3.18|) can be obtained to leading order in (3 2 from the Lagrangian 



- d^{x)d^{x) + ^ (cos p&{x) - 1). (3.21) 
2 p 



We argue that higher order corrections to the two-point Green function in a r (M 2 ) and 
to all orders in f3 2 should depend on the physical coupling constant a p h only 

-zA^\x,a r (M 2 )) = — ]l<Pyk(0\T(#(x)m--£i*(yk)'Mc = 

m ' J fc=i 

= -iA (m) (x,a ph ) (for m > 2). (3.22) 

In order to prove this assertion it is sufficient to analyse the renormalization of the causal 
two-point Green function to second order in a r (M 2 ) and to all orders in (3 2 . 

3.2 Two— point Green function to second order in a r (M 2 ) and to 
all orders in f3 2 

The correction to the two-point Green function to second order in a r (M 2 ) and to all 
orders in (3 2 is defined by 

-iA< 2 »(x,a r (M 2 )) = -I J J d 2 yi d 2 y 2 (o\T($(x)m ■ £int(yi) y- £int(y 2 ) m c = 

= _ - a 2 (M 2 ) V L_y_fl2(m-1) V- [ l ) o2(n 2 -l) 

2 -I ^(2.^ n ^ 2 (2n 2 )!^ 
x ff d 2 yi d 2 y 2 (O|T(i?(a;)0(O) : tf 2ni (l/i) :: ^(y 2 ) :)|0) c 



ni=2 V i; ' n 2 =l V 

x ff d 2 Vl d 2 y 2 (O|T(0(a;)0(O) : ^ 2ni (2/i) " :)|0> c 



2 rl A 1 J ^ (2ni)! P ^ (2n 2 )! P 

Airf 2 ^ (O|T(0(a;)0(O) : ^(2/1) " ^(2/2) :)|0) c . (3.23) 

This expression can be described in terms of two classes of topologically different Feynman 
diagrams depicted in Fig.l. In the momentum representation this correction is equal to 

-i A {2 \p,a r (M 2 )) =i \a r {M 2 ) Z x exp {\ f3 2 i A F (0, a r (M 2 ))' '" 



x 



00 

n=l 



2 

d 2 qi 1 /" d 2 q 2 



a r (M 2 )~p 2 



(2n) 2 i a r (M 2 ) - q 2 J (2vr) 2 i a r (M 2 ) - q 2 
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Figure 1: Feynman diagrams for corrections to the two-point Green function to second 
order in a and to arbitrary order in j3 2 . The left diagrams correspond to a non-vanishing 
expectation value for an even number of internal lines between the two vertices only, while 
the right diagrams describe a non-vanishing contribution for an odd number of internal 
lines between two vertices. 



x 



d 2 q 2n 1 1 

(2n)H a r (M 2 ) - q 2 n a r (M 2 ) - (p - q x - q 2 - . . . - q 2r , 



oo „ 
^4n + 2 / 

n=0 



d 2 qi 



d 2 q 2 



1 



x 



d 2 q 2 



n+l 



(2n) 2 i a r (M 2 ) - q\ J (2n) 2 i a r (M 2 ) 
1 1 



q\ 



(2vr) 2 2 a r (M 2 ) 
a r (M 2 ) 



qin+i 



a r (M 2 ) - (qi + q 2 + ■■■ + q2n+i) 



x 



x 



i \ a r (M 2 ) 



n=l 



M 2 
d 2 qi 



y 2 /8n- 


2 


r (-0 " 


2 






la r (M 2 ) -p 2 - 




1 




f d 2 q 2 


1 



(2tt) 2 2 a r (M 2 ) 
1 



ql 



(2?r) 2 ? a r (M 2 ) 
1 



ql" 



d 2 q 2 n 

(2n) 2 i a r (M 2 ) - q 2 n a r (M 2 ) - (p - q x - q 2 



q 2 r. 



oo „ 
n=0 J 



d 2 q 1 



d 2 q 2 



X 



d 2 q 2 



n+l 



(2n) 2 i a r (M 2 ) - q 2 J (2n)H a r (M 2 ) - q 2 " ' 
1 1 



(2tt) 2 z a r {M 2 ) - q 2 n+l a r {M 2 ) - (q x + q 2 + . . . + q 2n 



+1; 



(3.24) 



The common factor [a r (M 2 )Zi exp{j3 2 iAp(0, a r (M 2 )/2)}] 2 is caused by the summation 
of the infinite series of one-vertex-loop diagrams. Using (jl.4|) and ()3.19|) one can show 
that it is equal to a 2 h : 



a r (M 2 )Z 1 exp <^ - (3 2 iA F (0, a r (M 2 )) 



a r (M 



2j a r (M 2 )\P 2 /^ 2 



M 2 



a 



ph- 



Thus, we have shown that after renormalization the correction to the two-point Green 
function to second order in a r (M 2 ) and to all orders in the f3 2 is proportional to a 2 h . 
Then, since the two-point Green function, calculated to first in a r (M 2 ) and to all orders 
in j3 2 , is given by (|3.18jh in the correction of the second order a r (M 2 ) we can replace the 
two-point Green functions of the free sine-Gordon fields in the r.h.s. of ()3.24j) by (j3.18j) . 
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This gives 



-i A (2) (p,a r (M 2 )) =i«p h 



-i) 



x 



1 



d 2 q 1 



d 2 q 2 



(27r) 2 i a ph - q\ J (2n) 2 i a ph - q\ 



,2 ■ ■ - 



(2tt) 2 z a ph - g§„ a ph - (p - q 1 - q 2 



q2r 



CO „ 

+X> 4n+2 / 

n=0 J 



d'q! 



d 2 q 2 



(27r) 2 z a ph -qlJ {2ix) 2 i a ph - q 2 



2 • • • 



X 



<i 2 g 2 



n+l 



(27r) 2 z a ph - gf n+1 a ph - (gi + g 2 + • ■ • + g 2 



n+l, 



(3.25) 



In the momentum representation this proves relation ()3.22|) to second order in a r (M 2 ) 
and to all orders in (3 2 : 



l A^(p,a r (M 2 )) = - l A^ 2 \p,a ph ). 



(3.26) 



The proof of relation ()3.22|) to arbitrary orders in a r (M 2 ) and f3 2 demands only patience 
and perseverance. 



3.3 Non— trivial finite corrections to the dimensional coupling 
constant a p h 

As has been shown above quantum fluctuations around the trivial vacuum, calculated to 
first order in a r (M 2 ) and to all orders in /3 2 , lead to the renormalization of the dimensional 
coupling constant «o(A 2 ), which reduces to the replacement «o(A 2 ) — > a p h, where a p h is 
the physical (observable) dimensional coupling constant independent on both the cut- 
off A and the normalization scale M. In turn, quantum fluctuations around the trivial 
vacuum, calculated to second order in a r (M 2 ) and to all orders in /3 2 , induce non-trivial 
perturbative finite corrections to the physical coupling constant ct p h. 

The simplest correction of this kind is of order 0(a 2 h /3 2 ). It leads to the perturbative 
finite correction to the soliton mass, coinciding with that calculated by Dashen et al. 
017]. 



Keeping only the terms of order 0(f5 2 ) in (J3.25|) we get 



ia ph 



-i) 



la ph -p 2 - iO 

P 2 



d 2 qi 1 

{2n)H (a ph -q 2 -t0f 



An 



la ph -p 2 - iO 



(3.27) 



where we have taken into account that the integral over qi is equal to l/(47ra p h). 

The two-point Green function calculated to second order in a p h and to first order in 
f3 2 is equal to 



i A(p) 



-i 



« P h —jr — iU Air 



ttph — P i -I a e fr — p — iO' 



(3.28) 
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where a e g is defined by 

/ f3 2 \ 

«eff = « P h [1 - — J • (3.29) 

The effective Lagrangian of the sine-Gordon model defining the two-point Green function 
(13.28)) to leading order in (3 2 takes the form 

C(x) = - d^{x)d^{x) + ^ (cqbP#(x) - 1). (3.30) 

The soliton mass, calculated for the effective Lagrangian ()3.30|) . is equal to jnilZj(see also 
0) 

Ms - ~ ~' [ } 

This result coincides with that obtained by Dashen et al. [HJ E] ■ 

We would like to remind that the finite correction — y/a p h/^ has been interpreted in the 
literature as a singularity of the sine-Gordon model at [3 2 = 8n (see also 0). However, as 
has been conjectured by Zamolodchikov and Zamolodchikov [8,, such a singularity of the 
sine-Gordon model is superficial and depends on the regularization and renormalization 
procedure. This conjecture has been corroborated in [2]. 

In our present analysis of the sine-Gordon model the finite correction — ^/a^/n is 
a perturbative one. It is valid only for (3 2 87r and introduces no singularity to the 
sine-Gordon model at (3 2 = 8n. 



3.4 Non— trivial momentum dependent corrections to the two- 
point Green function 

Quantum fluctuations around the trivial vacuum, calculated to second order in a p h and 
to second order in (3 2 inclusively, lead to the non-trivial momentum dependence of the 
two-point Green function of the sine-Gordon field. From (|3.25|) we find the correction to 
the two-point Green function to order 0(a 2 h j3 4 ) inclusively. This reads 



d 2 qi 



i a 



" P h - V 
d 2 q 2 



d 2 Q! 

(2n) 2 i (a ph 
1 



1 



(27r) 2 z a ph -qfj {2n) 2 i a ph - q\ a ph - (p - q x - q 2 ) 



(3.32) 



The momentum integral of the contribution of order 0{f3 2 ) is equal to l/(47ra p h). For the 
calculation of the momentum integrals of the term of order 0(/3 4 ) we apply the Feynman 
parameterization technique. This gives 



d 2 q 1 



d 2 q 2 



(2tt) 2 z a ph -q 2 J (2vr) 2 i a ph - q\ a ph - (p - q x - 
dr)idr) 2 dr) 3 5(l -r} 1 -r] 2 - r] 3 ) 



Q2) 



1 r l 



167T 2 J J J a p h(^i^2 + V2V3 + Mi) + (~P 2 )%V2V3 
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1 



167T 2 
1 



drj / d£ 



x £n{ 



a P hV + (" P h - P 2 v) (! - V) £ C 1 - 

87r2 io ^("ph - P 2 ^)( 1 - ^)(4« P h^ + (« P h - P 2 ^)(l - ??)) 
/ ^{a pil -p 2 v){l + y/Aa ph ri + {a ph -p 2 r])(l^j) \ ^ 
W(«ph ~P 2 V){ 1 ~V) ~ V Aa phV + (« P h - p 2 v) (i - *7r 



We propose to analyse the behaviour of this integral in the asymptotic regime p 2 —>■ oo, 
where it can be calculated analytically. In this limit the main contribution to the integral 
over 7] comes from the domain rj ~ a v h/p 2 - Therefore, the integrand can be transcribed 
into the form 

d 2 q 1 1 f d 2 q 2 



(2n) 2 i a ph 



1 



1 



(27r) 2 z a ph 



1 



drj 



n 2 J (>ph - P 2 v) 



in 



p 2 T] 



ql a P h - (p 



a ph i] 



-£n 2 (^-) 
167r 2 p 2 Va p h^ 

i 



+ . . . . 



(3.34) 



The correction to the two-point Green function of order 0(a 2 h /3 4 ) inclusively, taken in 
the asymptotic regime p 2 — > oo, is equal to 



A {2 ' 2) (p,a ph ) 



i a 



P h 



r (-0 " 


2 (3 2 r 


-a P h -p 2 - 


Air . 



"ph 

47T p 2 



-ttph 



(3.35) 



Hence, the two-point Green function of the sine-Gordon field, accounting for the contri- 
butions of order 0(a 2 h f3 4 ) inclusively, takes the form 



A-\p) 



a 



P h 



47T 



^^£n 2 (V-) 
Vaph^ 



167T 2 p 2 



P 



(3.36) 



This shows that (i) all divergences can be removed by the renormalization of the dimen- 
sional coupling constant «o(A 2 ), (ii) the renormalized expressions are defined in terms of 
the physical coupling constant a p h, (iii) higher order corrections in a p h introduce a non- 
trivial momentum dependence and (iv) in the asymptotic limit p 2 — > oo the two-point 
Green function of the sine-Gordon field behaves as A(p) — > l/(— p 2 ). Such a behaviour 
is confirmed by the analysis of the two-point Green functions with the Callan-Symanzik 
equation (see Section 4). 



3.5 Physical renormalization of the sine— Gordon model 

Using the results obtained above we can formulate a procedure for the renormalization 
of the sine-Gordon model dealing with physical parameters only. Starting with the La- 
grangian (jl.lj) and making a renormalization at the normalization scale M 2 = a p h we 
deal with physical parameters only 

a ph = Z{\P 2 , a ph ; A 2 ) a (A 2 ), (3.37) 

where the renormalization constant Z\{fi 2 , a p ^; A 2 ) is equal to 

/ A 2 \ /3 2 /87r 

Z 1 (/5 2 ,« ph ;A 2 )= ( — ) . (3.38) 

VCKph/ 
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The renormalized Lagrangian is defined by 

C(x) = - {d^{x)d^{x) + ^ (cos^(x) - 1) + [Z x - 1) ^ (cos/3tf(a;) - 1) (3.39) 
2 p p 

with the renormalization constant given by Eq. (j3.37|) . From the relation ()3.20|) at M 2 = 
a p h one can obtain that 

a r (a P h) = a ph . (3.40) 

The calculation of perturbative corrections to the two-point Green function of the sine- 
Gordon model shows that the first order correction in a p h vanishes in accordance with 
Eq. ()3.17|) . Non-trivial corrections appear only to second and higher orders in a p h- 

4 Renormalization group analysis 

In this Section we discuss the renormalization group approach [H| [S] to the renormaliza- 
tion of the sine-Gordon model. We apply the Callan-Symanzik equation to the analysis 
of the Fourier transform of the two-point Green function of the sine-Gordon field. 

The Callan-Symanzik equation for the Fourier transform of the two-point Green func- 
tion of the sine-Gordon field (|3.5j) . which we denote below as —i A(p; a r (M 2 ), /3 2 ), is equal 
to 

" P- Yp + PMM 2 ),P 2 ) da ® {M2) - 2] A(p; a r (M 2 ),(3 2 ) = F(0, p; a r (M 2 ), (3 2 ), (4.1) 
where /3(a r (M 2 ), (3 2 ) is the Gell-Mann-Low function 

M^f± = f3MM 2 ),p 2 ). (4.2) 

The term 7(a r (M 2 ), (5 2 ) [Hj, describing an anomalous dimension of the sine-Gordon field, 
does not appear in the Callan-Symanzik equation (|4.1|) due to unrenormalizability of the 
sine-Gordon field d(x). The r.h.s. of (|4.1|) is defined by 

F(0,p; a r (M 2 ), (3 2 ) = jj d 2 xd 2 y e +i P ' x (0\T (e»(y)0(x)0(O) e % I ^ C ^(y)^j | ) c; (4.3) 
where £i n t{y) is equal to [2] 



Ant(y) = p(cos^(a;)-l). (4.4) 



Then, Q^(y) is the trace of the energy-momentum tensor Q IMU (x). For a (pseudo) scalar 
field $(x), described by the Lagrangian C(x), it is defined by [3| 

Using the Lagrange equation of motion one can show that 

^e^(z) = 0. (4.6) 
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For the sine-Gordon model the energy-momentum tensor M j,(x) reads 



e^(x) = d^{x)d v d{x) - g, u [^d x $(x)d x $(x) + p (cos/M(x) - 1) 

The trace of the energy-momentum tensor Mi ,(x) is equal to 

2a 



(3 2 



(cos/ft?(x) - 1) =2V[d(x)], 



(4.7) 



(4.t 



where V[i?(a;)] is the potential density functional of the sine-Gordon field i9(x). 

Since the trace of the energy-momentum tensor is proportional to the potential energy 
density, the Fourier transform F(0,p; a r (M 2 ),/3 2 ) can be related to the two-point Green 
function as 



F(0,p;a r (M 2 ),/r) = 2a r (M 



d 



da r (M 2 



-A(p;a r (M 2 ),(3 2 ), 



(4.9) 



where we have used the definition of the trace Q^(y) of the energy-momentum tensor 
Eq.flOD and the relation a = a r {M 2 ) Z x {f3 2 , M 2 ; A 2 ). 

Substituting (|4.9|) into ()4.1)1 we arrive at the Callan-Symanzik equation for the Fourier 
transform of the two-point Green function of the sine-Gordon field 



p 2 ^+(\p(a r (M 2 ),p 2 )-a r (M 2 : 



d 



da r (M 2 ) 



A(p 2 ;a r (M 2 ),(3 2 )=0, (4.10) 



where we have taken into account that A(p; a r (M 2 ), f3 2 ) should depend on p 2 due to 
Lorentz covariance. 

For the solution of (j4.10j) we have to determine the Gell-Mann-Low function (|4.2|) . 
For the coupling constant a r (M 2 ), defined by (J3.2(Jj) . the Gell-Mann-Low function is 



(3(a r (M 2 ),f3 2 ) = ^a r (M 2 ), 

47T 

where j3 2 = f3 2 /(l + /3 2 /8n) (|3.20|) . This gives the Callan-Symanzik equation 



(4.11) 



p 2 — 

dp 2 



+ 1 



0' 



•57T 



a T {M 2 







da r (M 2 



+ 1 



A(p> r (M 2 ),/^ 



(4.12) 



Setting A(p 2 ;a r (M 2 ),(3 2 ) = D(p 2 ;a r (M 2 ), (3 2 )/p 2 we get 



. ap A 



in 



a r {M 2 







da r (M 2 ) 



D(p 2 ;a r (M 2 ),(3 2 ) = 0. 



(4.13) 



Due to dimensional consideration the function D(p 2 ; a r (M 2 ) , (3 2 ) should be dimensionless, 
depending on the dimensionless variables p 2 = p 2 /M 2 and a = a r (M 2 )/M 2 , where M is 
a normalization scale. This gives 



p 2 — 



1 I a — — 



In J da 
14 



D(p 2 ;a } (3 2 ) = 0. 



(4.14) 



According to the general theory of partial differential equations of first order ^U|, the 
solution of (|4.14|) is an arbitrary function of the integration constant 



which is the solution of the characteristic differential equation 

^ dp 2 
p 2 



/ j3 2 \dp 2 da 
V 8irJ fj 2 a 



(4.15) 



(4.16) 



Hence, the Fourier transform of the two-point Green function of the sine-Gordon field is 
equal to 



h{p 2 ] a r {M 2 ) 1 p 2 ) = - 2 D 



a r (M 2 ) / p 2 \PV^- 
~~p~ 2 KM 2 ) 



(4.17) 



The argument of the D-function can be expressed in terms of the running coupling con- 
stant a r {p 2 ): 

^)= ar (M%^) =«*(-) (JL) =a pl ,(iL) . (4.18) 



The solution of the Callan-Symanzik equation for the Fourier transform of the two-point 
Green function of the sine-Gordon field is 

A(p ;a ph ,(3 ) = — D 



V 



(4.19) 



This proves that the total renormalized two-point Green function of the sine-Gordon field 
depends on the physical coupling constant a p h only. 

A perturbative calculation of the two-point Green function, carried out in Section 
3, gives the following expression for the function D[a r (p 2 ) / p 2 } in the asymptotic region 



oo: 



D 



-a r (p 2 ) - 


f a P h 




/? 4 


. p 2 


I p 2 


L 4tt 


16tt 2 



2 

' p N 



-1 



(4.20) 



Unfortunately, this is not able to reproduce the non-perturbative expression of the two- 
point Green function A(p 2 ; a p h,/3 2 ). 



5 Renormalization of Gaussian fluctuations around 
solitons 

We apply the renormalization procedure expounded above to the calculation of the con- 
tribution of quantum fluctuations around a soliton solution. We start with the partition 
function 



J Vd exp {i J d 2 x ^d^(x)d^(x) + p (cos/ft?(x) - 1) } = 

Wexpji J d 2 xC[d(x)}y (5.1) 
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Following Dashen et al. [HI El (see also [TT]) we treat the quantum fluctuations of the 
sine-Gordon field $(x) around the classical solution -&(x) = $ c i(x) + where (p(x) is 
the field fluctuating around ^(x), the single soliton solution of the classical equation of 
motion 



equal to El HH 



n #d(x) + j sin/M rf (x) = 



4 4 

— arctan(exp(- v /cio 7(x 1 - mi )) = — arctan(exp( v /cio a)), 



(5.2) 



(5.3) 



where u is the velocity of the soliton, a = 7(x 1 — mx°) and 7 = 1/yl — u 2 1 . 

Substituting #(x) = $ c e(x) + ip(x) into the exponent of the integrand of (|5.1|) and 
using the equation of motion for the soliton solution 1^(2) we get 



Z SG = exp \ % I d 2 x£[$ ce (x)} \ I Vip exp <j i / d 2 x [- d^{x)d^{x) 



a 



+ — sin/3$ rf (x) ((3<p(x) - sin/3yj(x)) + — cos (3$ cl >(x)(cos (3(p(x) - 1 



(5.4) 



Substituting (|5.3|) into (|5.4|) we obtain 



C[p(x)] 



a 



- d fl ip(x)d fM ip(x) + — (cos f3ip(x) - 1) 



2 

2a 1 
/3 2 cosh 2 ( v /ao 7(x 1 



[cosP<p(x) — 1) 

/3 2 cosh 2 ( v /ao 7(x x — mi )) ^ ^ ' ^ 
In terms of cx and r the exponent of the partition function (|5.4jl reads 



MX )) 

2c*o sinh( A /a^7(x 1 — mx )) 



(5.5) 



Z SG = exp ji y d 2 x£[$^(x)]} 



Of) 



TV ex P ^ * / drda - d^r, cr)9 M (^(r, cr) + — (cos/3^(r, a) - 1) 



L2 



2a r 



/3 2 cosh 2 (yjoi^a) 
where we have denoted 



cos 0<p(x) — 1) 



P 2 

f3 2 cosh 2 ( v /a^cr) 



(M^)-sinM^))]},(5.6) 



/ , , fd(p(r,a)Y fd(f(r, 



a. 



(5.7) 



1 In analogy with the "spatial" variable a we can define the "time" variable for the soliton moving with 
velocity u as t = j(x° — ux 1 ). In variables (t, cr) an infinitesimal element of the 2-dimensional volume 
d 2 x is equal to d 2 x = g?tcZct and the d'Alembert operator □ is defined by □ = d 2 /dr 2 — d 2 /da 2 . 
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The equation of motion for the fluctuating field (p(r, u) is equal to 

□y?(r,cx) + ^ sin (3p(r, a) = +2^ — ] sin pV(t,<t) 

p p cosh (y/aocr) 



a sinh(^cr) 
-2 — — 2 v (l-cos/fyfotr)). (5.8) 
p cosh (y^uj 

Dealing with Gaussian fluctuations only [HI Ej and keeping the squared terms in the 
Lagrangian ()5.5|) in the fluctuating field <p(x) expansion only, we transcribe the partition 
function 1)5.6)1 into the form 



Z S(: = exp<jz / d 2 x C{&ce(x)] 
1 



x / Vif exp 



% — J drda <p(r, a) 



□ + a 



2a 



cosh ( A /aocr) 



<p(r,a)}. (5.9) 



It is seen that ^Jo^Q has the distinct meaning of the mass of the quanta of the Klein-Gordon 
field <p(r, a) coupled to an external force described by a scalar potential 2 . 

Integrating over the fluctuating field (p(r, a) we transcribe the r.h.s. of (J5.9)) into the 
form 



Z SG = exp ji / d 2 x C[$ ce (x)] + z5<S[i^] }. 



(5.10) 



We have denoted 



exp{i5S[d ci ]} = exp{i J d 2 x 5C cS [d ce (x)]} 



Det(D + a ) 



exp | - ^ ^2 \ n + ^ y ^ y 



Det(n + a \ 

v cosh {^aoa) 



d 2 p 

W) 2 



£n(ao — p 2 



(5.11) 



where p is a 1+1-dimensional momentum. The second term in the exponent corresponds 
to the subtraction of the vacuum contribution. The effective action, caused by fluctuations 
around a soliton solution, is defined by 



5S[&) = J d 2 x5£ cS [$ C i(x)]=i^J2 inK + l J d * x J 



d 2 p 
{2ir)H 



in(a -p 2 ), (5.12) 



where A„ are the eigenvalues of the equation 

2a!o 
cosh 2 ( v /ao<j) 



(n + a 



) (p n (r, a) = A„v2„(r, a) 



(5.13) 



and (p n (r, a) are eigenfunctions. The quantum number n can be both discrete and contin- 
uous. This implies that the product over n in ()5.11)1 should contain both the summation 

2 The parameter ag should enter with the imaginary correction ag — > ccq — i 0. This is required by the 
convergence of the path integral 
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over the discrete values of the quantum number n and integration over the continuous 
ones. 

According to the Fourier method ^2] j the solution of equation (J5.13|) should be taken 
in the form 

<p n (T,a) = e - luJT ^ n (a), (5.14) 

where — oo < uj < + oo and i/j n (<j) is a complex function 3 . 
Substituting (pTTIJ) into lETTSI) we get 



d 2 +k 2 + 2a 



— )^n(r,a) = 0, (5.15) 



^da 2 cosh 2 ( v /a ^ 
where we have denoted 

k 2 = \ n + uj 2 - a . (5.16) 

This defines eigenvalues A n as functions of to and k 

X n = a - u 2 + k 2 . (5.17) 

The parameter k has the meaning of a spatial momentum — oo < k < +oo. The solutions 
of equation (|5.15|) are 4 



2 cosh( A /7i^cr) 



Ma) = i -ik + ^ tanh(^) | ^ %ka 
\J2tx y/k 2 + a 

where the eigenfunction ipb(cr) has eigenvalues X n = —uj 2 and the eigenfunctions ipk(v) 
have eigenvalues \ n = ao — uj 2 + k 2 . In the asymptotic region o — > oo the function ipk{&) 
behaves as 

1 e * + ^W, (5.19) 

V 27T 



where 8(k) is a phase shift defined by [TT] 

= 2arctan^P^. (5.20) 
The solutions ([5.18|) satisfy the completeness condition [TT] (see Appendix) 

dk 1>t(a 1 )ip k {<?) + M°')M°) = 5(a'-a). (5.21) 



3 Since ip n (T, a) is a real field, we have to take the real part of the solution (| B . 1 4|1 only, i.e. (p(r, a) — 
He (e ~ lWT ip((j)). Though without loss of generality one can also use complex eigenfunctions |71I3HTT]. 

4 The solutions of the equation Ij5.15|l are well-known ^TJ (see also [HE])- Nevertheless, we adduce 
the solution of this equation in the Appendix. 
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The fluctuating field (p(r, a) is equal to (see (A. 12)) 



<Pub(r, a) = —f= e ip b (a) = — = \ — e 

V27T V27T V 2 cosh( v /a^c7) 



V«fc(T,ffJ = -j==e Wk(v) = - — e . (5.22) 



In terms of the eigenvalues X n = — uj 2 and A n = a — cu 2 + /c 2 and eigenf unctions ()5.22|) 
the effective action 5S is determined by 

i r f+oo j f+oc 

5S[$ ce ] = --J d 2 x J —J dk\^ k (x)\ 2 £n(a -uj 2 + k 2 ) 
-- [ d 2 x — \4, b ( x )\ 2 ln{-u 2 ) + - [ d 2 x [ -fl-£n(a - P 2 ). (5.23) 



2m Ir K n v ' 2 J J (2yr) 2 i 

Using the explicit expressions for the eigenfunctions ipk{x) and ipb{x) we reduce the r.h.s. 
of (ED) to the form 

do; f +OQ dk k 2 + a tanh 2 ( v / aocr) 2 2 



I /" ^ T 00 ^ V^£_^ ^(^2^1 /" ^ /" J^P £n ( ao _ p 2y (5 24) 



2ttz 2 cosh 2 (^cr) v ' 2 J J (2tt) 2 « 
Introducing the notation 

r^T= ; = 5 (1 - <x»P0 u {x)) = f- V[d ct {x)] (5.25) 
cosh (i/aoO") ^ 2ao 

we obtain the effective Lagrangian 8C e s\&d{ x )]- It is equal to 5 

1 r +co dij) f +oc dk 1 

5£ cS [Mx)] = --P 2 V[$ c e(x)] ^ — yy ; [M~^) - M«o - ^ 2 + fc 2 )] 

4 J_ 00 2m 2tt k 2 + a 

1 /"+ 00 du [ +co dk n , , , 2N 1 /" d 2 p n . 2N , 

' ' £n{a -uj 2 + k 2 ) + - —^-£n{a -p 2 ). (5.26) 



2y_ 00 27riy_ 0O 2vr v y 27 (2vr) 2 z 

The two last terms cancel each other. This gives 

1 f +oc dt,i f +oc dk 1 

6£ eS [# ce (x)} = —tfV^ix)] / / -^^[£n(-uj 2 ) - £n(a - uj 2 + k 2 )]. 

(5.27) 

After the integration by parts over uj the effective Lagrangian <5£ e ff [^(^O] is defined by 

1 /" + °° f +oc din 1 
<5£ eff [^(x)] = -/5V[^(x)] / ^/ 2 ■ ?2 - ■ (5-28) 



2 L v y V-oo 27ry_ ao 2vria -w 2 + A; 2 -iO' 



5 In the contribution of the mode A n = — lj 2 we have used the integral representation 

' 50 dfc 1 



2\/<^o V-oo 27r fc 2 + a 
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The appearance of the imaginary correction —% is caused by the convergence of the path 
integral Q PS- 

Integrating over u we reduce the r.h.s. of (J5.28j) to the form 



dk 



4tt Veto + k 2 
(3 2 Trrn /- +00 dA; y : , 9 ^(fc) 

dk 



/ dk 
— V^Ttf 
-oo 47r 



(5.29) 



where 8(k) is a phase shift defined in Eq. ()5.2()|) . We discuss this expression in Section 7 
in connection with the correction to the soliton mass caused by Gaussian fluctuations. 

For the analysis of the renormalizability of the sine-Gordon model the momentum 
integral in the effective Lagrangian 5£ e g {& c e{x)} should be taken in a Lorentz covariant 
form (J5.28J) and regularized in a covariant way. Making a Wick rotation uj — ■> iuo and 
passing to Euclidean momentum space we define the integral over u and k in (|5.28|) as [2 



dk 
2t7 



duo 



2m «o 



uj- 



1 = ^tn(^\ 

+ k 2 -i0 4vr \a J' 



(5.30) 



where A is an Euclidean cut-off j2]. The effective Lagrangian, induced by Gaussian 
fluctuations around a soliton solution, is equal to 



a 



5£ cS [$ C i(x)] = — 



A 2 M 



57T 



(cos/3^(x) - 1). 



(5.31) 



The Lagrangian ()5.31|) has the distinct form of the correction, caused by quantum fluc- 
tuations around the trivial vacuum calculated to first orders in a (A 2 ) and (3 2 . 

The total Lagrangian, accounting for Gaussian fluctuations around the soliton solution 
amounts to 



1 /32 A 2 

C eS (x) = -d^ ce (x)d^ ce (x) + p [l - |-£n(— ) 



(cos/3^(x) - 1). (5.32) 



This coincides with Eq.(6.7) of Ref.|2]. 

As has been shown in the dependence of the effective Lagrangian ()5.32j) on the 
cut-off A can be removed by renormalization with the renormalization constant (jl.4|) 



o 



8-7T V a ' 



a r {M 2 



a r (M 2 )Z 1 
A 2 \ 



8tt \a r {M 2 )Z l 
A 2 



2 

—in 
8tt \a r (M 2 



a r (M 2 



+ 8tt \ M 2 



, (5.33) 



where we have kept terms of order 0(/3 2 ) in the /3 2 -expansion of the renormalization 
constant (jl.4|) . This gives the effective Lagrangian 



-d^ ce (x)d^ ce (x) + 



a T {M 2 
~1 2 ~ 



1 + t in f^p 
8vr V M 2 



(cos/?^(x) - 1). (5.34) 
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We can replace the coupling constant a r (M 2 ) by the physical coupling constant a p h, 
related to a r (M 2 ) by (l3~T31) 



a 



P h 



a r {M 2 



+ 8tt I M2 



(5.35) 



where we have kept terms of order 0(/3 2 ) only. Substituting ()5.35j) in ()5.34j) we get 

1). (5.36) 



4ffV) = \d^{x)d^ ct {x) + ^ (cos/M c ,(x) 



The effective Lagrangian (|5.36|) coincides with the Lagrangian, renormalized by the quan- 
tum fluctuations around the trivial vacuum (|3.21|) . and corroborates the result obtained 
in 0. 

We would like to emphasize that analysing the renormalization of the sine-Gordon 
model, caused by Gaussian fluctuations around a soliton, one can see that Gaussian 
fluctuations are perturbative fluctuations of order O{a r {M 2 )0 2 ) valid for 1 -C 8tt, which 
cannot be responsible for non-perturbative contributions to the soliton mass at (3 2 = 8ir. 



6 Renormalization of the soliton mass by Gaussian 
fluctuations in continuous space— time 

Using the effective Lagrangian Eq. (j5.29|) one can calculate the soliton mass corrected by 
quantum fluctuations. It reads 



P 2 

where AM S is defined by 

AM S = - £°° dx 1 SC^eix 1 )] = J +C ° g A + (6-2) 

This corresponds to the correction to the soliton mass, induced by Gaussian fluctuations, 
without a surface term —^/oq/tt jT4"| IT5]. 

In the Lorentz covariant form the correction to the soliton mass reads 



/+oo 
dx 1 SCrtficiix 1 )] 
-co 



+°° duj f +oc dk r* / 2^ „ , 2 , 2m 

/ — , n v \£n(-uj 2 ) - £n(a - uo 2 + k 2 )} 

^ 27riJ_ 00 2nk 2 + a l 1 ; 1 n 

^_ f +co dk f +co du 1 . . 

= -2^ / — / — 6.3 

y_oo 2vr > /_ 00 2m a - uj 2 + k 2 - % 

where we have taken the effective Lagrangian defined by ()5.27j) and integrated over u 
by parts. Using the result of the calculation of the integral ()5.30|) we get the following 
expression for the soliton mass corrected by Gaussian fluctuations 



_ 8y/a^) ygggh j A 2 
s ~ P 2 2vr la (A 2 ) 
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(6.4) 



The removal of the dependence on the cut-off A should be carried out within the renor- 
malization procedure. 

Replacing ao(A 2 ) in (J6.4)) by a r (M 2 )Zi(f3 2 , M 2 ; A 2 ), where the renormalization con- 
stant Z x {f3 2 , M 2 ; A 2 ) is defined by Eo. lfOj) . we get 



M _ 8^a r (M 2 )Z 1 ^a r {M 2 )Z x ^ 



ff 1 



2tt 



A 2 



a r (M 2 )Z 1 



(6.5) 



The renormalization constant Z\ should be expanded in power of f3 2 to order 0(/3 2 ). This 
gives 



B 2 ( A 2 \ 
1 8tt \M 2 J 

Substituting ()6.6|) into ()6.5|) and keeping only the leading terms in 8 2 we get 
M,= ^,/*OT l + £fa(^ 



Using Eq. (j5.35|) we can rewrite the r.h.s. of ()6.7j) in terms of a^: 



/3 2 



(6.6) 



(6.7) 



(6.8) 



The mass of a soliton M s depends on the physical coupling constant a p h- Hence, the 
contribution of Gaussian fluctuations around a soliton solution is absorbed by the renor- 
malized coupling constant a p h and no singularities of the sine-Gordon model appear at 
(3 2 = 8tt. 

This result confirms the assertion by Zamolodchikov and Zamolodchikov jS], that the 
singularity of the sine-Gordon model induced by the finite correction — ^/a^h/V to the 
soliton mass, caused by Gaussian fluctuations around a soliton solution, is completely due 
to the regularization and renormalization procedure. This has been corroborated in j2|. 

As has been shown in Section 3, a non-trivial finite correction to the soliton mass 
appears due to non-Gaussian quantum fluctuations of order of ctp h /3 2 (see Eq. (j3.31|) ): 



M, 



~T 2 ~' 



(3.31) 



The second term in Eq. (j3.31j) coincides with that obtained by Dashen et al. |EIIZ|- How- 
ever, it is meaningful only as a perturbative correction for /3 2 -C 87r, which cannot produce 
any non-perturbative singularity of the sine-Gordon model at f3 2 = 8tt. 

We have obtained that the soliton mass M s does not depend on the normalization 
scale M. This testifies that the soliton mass M s is an observable quantity. 



7 Renormalization of soliton mass by Gaussian fluc- 
tuations. Space— time discretization technique 

Usually the correction to the soliton mass is investigated in the literature by a discretiza- 
tion procedure [Hj (see also [IS])- The soliton with Gaussian fluctuations is embedded 
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into a spatial box with a finite volume L and various boundary conditions for Gaussian 
fluctuations at x — ±L/2. In such a discretization approach time is also discrete with a 
period T, which finally has to be taken in the limit T —>■ oo. The frequency spectrum 
is u m = 2nm/T with m = 0, ±1, . . .. For various boundary conditions spectra of the 
momenta of Gaussian fluctuations around a soliton and of Klein-Gordon quanta, cor- 
responding to vacuum fluctuations, are adduced in Table 1. According to Table 1, for 
various boundary conditions the corrections to the soliton mass are given by 



AMj p) 

oo 

E K 



lim lim 

T^oo L^oo 2iT 



oo oo 
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£n(a 



4ir 2 m 2 
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+ A 2 (L) -£n[a 



A-n 2 m 2 

J^2 



+ *)]}, 



(7.1) 



m=— oo 



a ° ^ at L — > oo, and the abbreviations (p), (ap) and (rw) mean 



where A 2 (L) ~ a$e 
periodic, anti-periodic boundary conditions and rigid walls, respectively. 

For the summation over m we use the formula, derived by Dolan and Jackiw |16j : 
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m=— oo 



\ln 



4m 2 7r 2 

_ r p2 



+ a< 



£n( 



4m 2 7r 2 

^ _J^2 



i(a-b)T + 2£n 
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-ibT 



(7.2) 



Taking the limit T — ► oo we arrive at the expression 

oo 

n=l 

oo 

n=2 
^ oo 



AM, 



h hm < 

2 L^oo 



periodic BC 

anti — periodic BC (7.3) 
rigid walls, 



n=2 



where BC is the abbreviation of boundary conditions. 

The aim of our analysis of AM S within the discretization procedure is to show that 
the discretization procedure gives AM S in the form of ()6.3j) . 

The subsequent analysis of AM, we carry out for periodic boundary conditions only. 
One can show that for anti-periodic boundary conditions and rigid walls the result is the 
same. 
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Table 1: The spectra of the momenta of Gaussian fluctuations around a soliton and the 
Klein-Gordon quanta. The modes, denoted by (*) are due to the bound state. 



PERIODIC BC: 



Soliton Sector 



Vacuum Sector 



k n L + S(k n ) = 2nn 
n = : (*) 

n = 1 : qi = £ + 0(L" 



lx 
2x 
2x 



q n L = 2rm 
n = : qo = 
n = 1 : q ± = ^ 



Eoo re 
n=l ~ 



+0(L- 2 ) 2vr 



rf<5(fc) \ 



Eoo roo d(J r 

n=l ~^ J^L 



ANTI-PERIODIC BC: 



Soliton Sector 



Vacuum Sector 



k n L + 5(k n ) = (2n - 1)tt 
n = 1 : (*) + ki = 



(l + l)x 
<- 2x -> 



g n L = (2n - l)7r 
n = 1 : gi = f 



OO . poo 

71=2 



+0(L- 2 ) 2tt 



dk J 



Eoo 
n=2 



2^ 



RIGID WALLS: 



Soliton Sector 



Vacuum Sector 



k n L + 5(k n ) = nn 
n — 1 : (*) 



lx 
lx 



q n L = nn 

n = \:q x = l 



Y^oo roo dk f t , dS(k) \ 

l^n=2 ~* J2L +0 (L- 2 ) n \ ^~ dk ) 



Eoo 
n=2 



J -r- TV 



24 



For the next transformation of the r.h.s. of (|7.3|) we propose to use the following 
integral representation 



a + k 2 n - u? 



«o + q. 



2 -LU 2 



This gives 
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+00 
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(7.4) 



(7.5) 



To the regularization of the sum over n we apply the mode-counting regularization pro- 
cedure [T3] : 
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(7.6) 



Passing to the continuous momentum representation (T^j and using the spectra of the 
momenta of Gaussian fluctuations and the Klein-Gordon fluctuations (the vacuum fluc- 
tuations) adduced in Table 1 we transcribe the r.h.s of ()7.6|) into the form 
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where the limits are equal to (see Table 1): 
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Rearranging the limits of integrations we get 



AM„ (p) 



+ 



lim lim 

^ 27T* L^oo N->oo 



n/L 



dk dS(k) 
2tt dk 



£n(a + k 2 - J 2 - iO) 



+ 



lim L 

27TZ L^oo 



00 
+00 



^-in{a^ + k 2 - uj 2 - iO) 

n/L 271 



— : lim lim L 

2ni L-^oo N-^oo 



in 



fejv 



2^ 



£n(a + k 2 - cj 2 - z'O) 



25 



2 



hm / f-t-Mao + k -aT-zO) 

2m l^oo J n , L 2% dk 



.2 , .2 



+ / lim L / — £n(a + k - u - i 0) 

2m l-+oo J w/L 2vr 1 ' 

_ /" +00 ^_ i im i im L kN ~ qN in{ ao + k* 2 -uj 2 -zO). (7.9) 

For the last term we have applied the mean value theorem with — 5(qi\r)/L < k* < qjy. 
Since the difference k^ — qN = &{q.n)/L = ^/cTq/t^N is of order 0(1/N), this term vanishes 
in the limit iV — > oo 6 . As a result we get 

AM p) = -1—1 + / 1 m / ^ £n(a + k — u> — iO) 

+ / ■— lim L / — M«o + fc 2 - cu 2 - z 0). (7.10) 



2717 £— >oo y^/x, 27T 

Taking the limit L — > oo and applying to the last term the mean value theorem we arrive 
at the expression 



, x duj dk dS(k) „ , ,,, .„. ,„,,, 



which we transcribe into the form 



a n jr(ri\ v/«o r°° duo f°° dk d8(k) n . 2 . n s f + °° du n , 2 . n . 



hOO 



+ ^jf l^^ + t'-^-iOJ-^-iO)]. (7,2) 

The next steps of the reduction of the r.h.s. of (J7.12j) to the form (J6.3|) are rather straight- 
forward. First, one can easily show that 

+°° d00 f°° dk d8(k) „ , 2 f + °° du „ , 2 

= /" +0 ° [£ n (« _ ^ _ i 0) _ £n( _^2 _ • Q)] = v^o (7 13) 

J-oo 4vrz 2 

and, second, integrating over uj by parts the last integral in ()7.12|) can be reduced to the 
form 

AMi p ) = -2^ / — / — 2 ^, 2 .. . 7.14 

J_ 00 2ttJ_ oc 2ma -uj 2 + k 2 -i0 



6 We would like to emphasize that exactly the term of this kind leads to the contribution of the finite 
surface term — ^fag/ir in a regularization procedure using the expressions l|7.3|) without turning to the 
integral representation l|7.4|l . 
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Thus, we have shown that the correction to the soliton mass, induced by Gaussian fluc- 
tuations around a soliton and calculated by means of the discretization procedure, agrees 
fully with that we have obtained in continuous space-time (|6.3|) . 

Hence, the renormalization of the soliton mass, caused by Gaussian fluctuations calcu- 
lated within the space-time discretization technique, coincides with the renormalization of 
the soliton mass in continuous space-time. We would like to emphasize that the obtained 
result 



M4j) does not depend on the boundary conditions. 
The calculation of the functional determinant within the discretization procedure has 
confirmed the absence of the correction — a/oo/ 7T. This agrees with our assertion that such 
a correction does not appear due to Gaussian fluctuations around a soliton, corresponding 
to quantum fluctuations to first orders in ao(A 2 ) and (3 2 . 

The reduction of AM s (p) of Eq-fTTJ) to the expression ()6.3|) can be carried out directly. 
First, summing over n within the mode-counting regularization procedure and taking the 
limit L — > oo we arrive from AM^ of Eq. fj7.lj) at 
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(7.15) 
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Now we use the integral representation 5 and get 
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s T^ociT ^ J 2n a + k 2 L V T 2 / V T 2 )\ 
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For the other boundary conditions we get the same result. 

Thus, we have shown that the discretized version of the correction to the soliton mass 
reduces to the continuum result if one transcribes first the sum over the quantum number 
n of the momenta of Gaussian and vacuum fluctuations into the corresponding integral 
over the momentum k. 



8 Conclusion 

We have investigated the renormalizability of the sine-Gordon model. We have analysed 
the renormalizability of the two-point Green function to second order in a and to all 
orders in (3 2 . We have shown that the divergences appearing in the sine-Gordon model 
can be removed by the renormalization of the dimensional coupling constant cto(A 2 ). We 
remind that the coupling constant j3 2 is not renormalizable. This agrees well with a 
possible interpretation of the coupling constant [3 2 as h [TJ H2j. The quantum fluctua- 
tions calculated to first order in a r (M 2 ), the renormalized dimensional coupling constant 
depending on the normalization scale M, and to arbitrary order in f3 2 after removal of 
divergences form a physical coupling constant a p h, which is finite and does not depend on 
the normalization scale M. We have argued that the total renormalized two-point Green 
function depends on the physical coupling constant a p h only. In order to illustrate this 
assertion (i) we have calculated the correction to the two-point Green function to second 
order in a r (M 2 ) and to all orders in f3 2 and (ii) we have solved the Callan-Symanzik 
equation for the two-point Green function with the Gell-Mann-Low function, defined to 
all orders in a r (M 2 ) and f3 2 . We have found that the two-point Green function of the 
sine-Gordon field depends on the running coupling constant a r (p 2 ) = « p h(p 2 /tt P h)^ 2//87r , 
where f3 2 = (3 2 /{l + (3 2 /8n) < 1 for any (3 2 . 

We have shown that the finite contribution of the quantum fluctuations, calculated to 
second order in a p h and to first order in f3 2 , leads to an effective Lagrangian 

C(x) = \ d^{x)d^{x) + ^ (cos/M(x) - 1) (3.30) 

2 p 

with the effective coupling constant a e s = a p h (1 — /? 2 /4tt) , which defines the soliton mass 
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The term — ^Japh/ V is a perturbative finite correction to the soliton mass, which coincides 
with the result obtained by Dashen et al. [HI El • Such a correction has been interpreted 
in the literature as a non-perturbative contribution leading to a singularity of the sine- 
Gordon model at [3 2 = 8tc (see also Pj). However, in our approach such a correction is 
a perturbative one valid for f3 2 <C 8n and does not provide any singularity for the sine- 
Gordon model in the non-perturbative regime f3 2 ~ 8tt. This confirms the conjecture 
by Zamolodchikov and Zamolodchikov [Hj that a singularity of the sine-Gordon model at 
/3 2 = 87r [HllZlini is a superficial one and depends on the regularization and renormalization 
procedure. This has been corroborated in |2j. 

In addition to the analysis of the renormalizability of the sine-Gordon model with 
respect to quantum fluctuations relative to the trivial vacuum, we have analysed the 
renormalizability of the sine-Gordon model with respect to quantum fluctuations around 
a soliton. Following Dashen et al. [HI El and Faddeev and Korepin jU] we have taken into 
account only Gaussian fluctuations. 

For the calculation of the effective Lagrangian, induced by Gaussian fluctuations, 
we have used the path-integral approach and integrated over the field ip(x), fluctuat- 
ing around a soliton. This has allowed to express the effective Lagrangian in terms of 
the functional determinant. For the calculation of the contribution of the functional de- 
terminant we have used the eigenfunctions and eigenvalues of the differential operator, 
describing the evolution of the field <p(x). We have shown that the renormalized effective 
Lagrangian, induced by Gaussian fluctuations around a soliton, coincides completely with 
the renormalized Lagrangian of the sine-Gordon model, caused by quantum fluctuations 
around the trivial vacuum to first order in and to second order in f3 2 . After the removal 
of divergences the total effective Lagrangian, caused by Gaussian fluctuations around a 
soliton, is equal to (|3.21|) 

C{x) = - d^(x)d^{x) + ^ (cos/M(x) - 1). (3.21) 

2 p 

This implies that Gaussian fluctuations around a soliton do not produce any quantum 
corrections to the soliton mass. After the removal divergences the soliton mass is equal 
to the mass of a soliton, calculated without quantum corrections, with the replacement 
«o —* a ph- Hence, no non-perturbative singularities of the sine Gordon model at (3 2 = 8n 
can be induced by Gaussian fluctuations around a soliton. 

For the confirmation of our results, obtained in continuous space-time, we have calcu- 
lated the functional determinant caused by Gaussian fluctuations around a soliton within 
the discretization procedure with periodic and anti-periodic boundary conditions and 
rigid walls. We have shown that the result of the calculation of the functional determi- 
nant (i) coincides with that we have obtained in continuous space-time and (ii) does not 
depend on the boundary conditions. 
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Appendix. Solutions of the differential equation for 
Gaussian fluctuations around a soliton 



Making a change of variables £ = tanh( v /ao<j) we reduce (J5.15|) to the form 



d 



#(0 

di j 



+ 



s(s + l) 



m = o, 



(a.i; 



where s(s + 1) = 2 and e 2 = —k 2 /a = 1 — oo 2 /a (Hi- 
Substituting = (1 — £ 2 ) e ' 2 w(£) and denoting u = (1 — £)/2 we arrive at the 
equation 



u(l - u)w" + (e + 1)(1 - 2u) w' - (e - s)(e + s + 1) w = 0. 



(A.2) 



The solution of this equation can be given in terms of hypergeometric functions F(a, b; c; z) 
HS1 (see also [H]). It reads 



w(0 = + = d F(e - s, e + s + 1; e + 1; LA) 



+C 2 



1-^ 



s, s + 1; 1 - e; 



(A.3) 



where Ci and C 2 are the integration constants. 

The parameter s acquires two values s = —2, +1, which are solutions of the equation 
s(s + 1) = 2. Since for both case the hypergeometric functions coincide [TH], setting s = 1 
we obtain 



w(£) = wW(£) + w^{0 = d F[ e - 1, e + 2; e + 1 



,(2) 



"Co 



F[ -l,2;l-e: 
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For arbitrary e the solution u/ 2 )(£) can be reduced to the polynomial ^ 
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^ {2) (0 = c 2 



F( -l,2;l-e: 



1-e 



Co 
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(A.5) 



One can show that the solution w^'(^) can be reduced to a polynomial too. For this aim 
we suggest to use the following relation for hypergeometric functions 



F[a,b: c; 



1 -£\ _ T(c)T(c-a-b) 
~ T(c-a)T(c-b) 



F[a,b;a + b — c+l; 
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i + Ac-o-6r(c)r(o + 6-c) 



r(a)r(6) 



F(c — a, c — 6:c— a — fe+l; 



2 



(A.6) 



For the solution w^($) we get 
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^ (0 = CiF^e- l,e + 2; e + 1; = Ci — ±gg= £ F ( e - 1, e + 2; e + 1 ; l±i 

r(£ + i)r(e) + / i + e\ _ r /i+£\-e + £ 

+Cl r( e -i)r( e + 2) l—J ^l 2 '" 1 ' 1 " 6 '^-;-^^; 1+1' (A - 7) 

Thus, the fluctuating field ip(£) is equal to 

✓ _£ve/2 e _|_£ /Uf v /2f_ £ 

*tt> = °> 2 "(i+f) T^ 2 '(t^) 137- < A ' 8 > 

In terms of a the classical solution for the fluctuating field reads 



^j(a) =d(e + tanh( Vcwr)) e + C 2 (-e + tanh( v ^a)) e + e v / ^ CT , (A.9) 



where we have redefined the integration constants. These solutions describe a bound state 
for e = ±1 and a scattering state for e = ±i/c/y / Oo. They are 



2 cosh( v /«^(T 



^(a) = 4= — + ^ tanh ^) e (A.10) 
V27T y/k 2 + a 

The wave functions V , fc( cr )j given by (A. 10), are normalized to the ^-function as 

30 1 / d d \ +°° 

daip* k ,{a)il; k (a) = — _ fc2 ^*,(<r) — ^(cr) - M*)^ i ) t'( cr )) ^ = ~ *0- 

(A.ll) 

The fluctuating field <^(r, a) is equal to 



<Pwb(r, a) = —== \ ^- — e 

y/2n V 2 cosh^^/a^cr) 



/ \ i ~ ik + v^o tanh( v ^q) + ^ a 

y W fcCr,o-; = - e . (A.12) 

They are normalized by 

OO f+OO 

dri Pl'bi. T ^) l P^b{.r,o) = S(u' - lu), 



-OO J — OO 
-OO f + OO 



drda <fil,r k >(r, cr)^ wfc (r, cr) = - d>)5(A;' - fc). (A.13) 



OO »/ — OO 



The eigenfunction ipb{°') h as the eigenvalue zero, c<j = 0. This corresponds to the bound 
state [TT]. The eigenfunction ^(cx) is normalized to unity 

Ax|<Ma)| 2 = l. (A.14) 

) 
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The solutions (A. 10) satisfy the completeness condition [TT] 



+00 



dk r k (° ' )M°) = S(a'-a). (A.15) 



The proof of the completeness condition (A.15) 

^+00 r +oo dJ . 



2tt 



[tanh(v^a) -tanh( v ^a')] / ^ e ik ( a a ') 

+a [tanh( v ^a')tanh( v ^a) - 1] / tr TT^ e lk (° ~ °') 

2ir k z + a 



I 1 1 (A . 16) 

2 cosh( v /aocr / ) cosh( v /aocr) 



Integrating over we get 



dk r k (o- ' )^(a) + V 6 (o- ' )^(*) = ' - *) V " " 



2 cosh( v /a T o<j / ) cosh( v /aoa) 
x e - v^ok " CT ' I [ e (o- _ ^ ' ) sinh(^(<7 - a ' )) + cosh(vM CT - a ' ))] 



+ ^ -. j— . - j— . = £(a'-<r), (A.17) 
2 cosh( v /aocr ) cosh^y^a) 



where e(cr — cr ' ) is the sign-function. The term, proportional to e ~v a ^\ a ~~ a I, is given 
by the contributions of the second and the third terms of the r.h.s. of (A. 16). For the 
derivation of (A.17) we have used the relation 

e(a-a') smh(^((a - a ' )) + cosh( v ^(a - a ' )) = e ~ ' I . (A.17) 



Thus, the contribution of the second and the third terms in the r.h.s. of (A. 16) cancel 
the contribution of the zero-mode ipb{c)- This completes the proof of the completeness 
condition (A.15). 



32 



References 

[I] S. Coleman, Phys. Rev. D 11, 2088 (1975). 

[2] M. Faber and A. N. Ivanov, J. Phys. A 36, 7839 (2003) and references therein. 

[3] C. Itzykson and J.-B. Zuber, in QUANTUM FIELD THEORY, McGraw-Hill, New 
York, 1980. 

[4] M. E. Peskin and D. V. Schroeder, in QUANTUM FIELD THEORY, Perseus Books, 
Reading, Massachusetts, 1995. 

[5] S. Weinberg, in THE QUANTUM THEORY OF FIELDS, Vol. II Modern Applica- 
tions, Cambridge University Press, 1996. 

[6] R. F. Dashen, B. Hasslacher, and A. Neveu, Phys. Rev. D 10, 4130 (1974). 

[7] R. F. Dashen, B. Hasslacher, and A. Neveu, Phys. Rev. D 11, 3424 (1975). 

[8] Alexander B. Zamolodchikov and Alexey B. Zamolodchikov, Ann. of Phys. 120, 253 
(1979). 

[9] L. D. Faddeev and V. E. Korepin, Phys. Rep. 42, 1 (1978). 

[10] R. Courant and D. Hilbert, in METHODS OF MATHEMATICAL PHYSICS, Partial 
Differential Equations, Vol. II, Interscience Publishers, John Wiley & Sons, New York, 
1962. 

[II] J. Rubinstein, J. of Math. Phys. 11, 258 (1970). 

[12] N. N. Lebedev, I. P. Slalskaya, and Y. S. Uflyand, in Problems of Mathematical 
Physics, Prentice-Hall, Inc., Englewood Cliffs, New Jersey, 1965, pp. 55-102. 

[13] (see jU p.286) 

[14] R. Rajaraman, in Solitons and Instantons, North-Holland, 1982. 

[15] A. Rebhan and P. van Nieuwenhuizen, Nucl. Phys. B 508, 449 (1997). 

[16] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974) 

[17] M. Faber and A. N. Ivanov, Eur. Phys. J. C 20, 723 (2001); hep-th/0105057 

[18] L. D. Landau and E. M. Lifshitz, in QUANTUM MECHANICS, Pergamon Press, 
New York, 1965. 

[19] N. N. Lebedew, in SPEZIELLE FUNKTIONEN UND IHRE ANWENDUNG, Wis- 
senschaftsverlag, Bibliographisches Institut, Mannheim, 1973. 

[20] HANDBOOK OF MATHEMATICAL FUNCTIONS, with Formulas, Graphs, and 
Mathematical Tables, ed. by M. Abramowitz and I. E. Stegun, U.S. Department of 
Commerce, National Bureau of Standards, Applied Mathematics Series • 55, 1972. 



33 



